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Starting from a classical picture of shear viscosity we construct a stationary veloc- 

(N 

ity gradient in a microscopic parton cascade. Employing the Navier-Stokes ansatz we 
extract the shear viscosity coefficient r]. For elastic isotropic scatterings we find an 
■ excellent agreement with the analytic values. This confirms the applicability of this 

method. Furthermore for both elastic and inelastic scatterings with pQCD based 
cross sections we extract the shear viscosity coefficient r\ for a pure gluonic system 
and find a good agreement with already published calculations. 
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PACS numbers: 47.75. +f, 12.38.Mh, 25.75.-q, 66.20.-d 



INTRODUCTION 



Recent results of the Relativistic Heavy Ion Collider (RHIC) and of the Large Hadron 
Collider (LHC) indicate the formation of a new state of matter, the quark-gluon plasma 



(QGP), in relativistic heavy-ion co 
v 2 observed in these experiments 



lisions. The large value of the elliptic flow coefficient 



Hi 



a nearly perfect fluid. This has been confirmed. 
[5 11 1 and microscopic transport calculations 



44| leads to the indication that the QGP behaves like 
jy calculations of viscous hydrodynamics 
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131 ] . However, the shear viscosity coeffi- 



cient r] has a finite value, possibly close to the conjectured lower bound rj/s = 1/4-7T from 
the correspondence between conformal field theory and string theory in an Anti-de-Sitter 
space jlj]. In comparison to ideal bydrodynamic calculations [15], dissipative hydrodynamic 



formalisms with finite rj/s ratio [5 - 1 1 1 demonstrate a better agreement of the differential 



elliptic flow V2(pt) with experimental data. The shear viscosity is therefore an important 
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parameter in viscous hydrodynamics but needs to be calculated from microscopic theory 
The rj/s ratio was obtained in a full leading order pertubative QCD calculation in 



Ref. |16j . The Boltzmann-Vlasov equation and quasi-particle picture were recently em- 
ployed to calculate the 77/s ratio of a gluon gas in Ref. 12| • The shear viscosity coefficient 
has also been extracted from microscopic transport calculations with BAMPS (Boltzmann 
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19| using expressions based on a 
20H and the entropy principle 



Approach of Multi Parton Scatterings) simulations 
first-order gradient expansion of the Boltzmann Equation 
underlying the second-order Israel- Stewart hydrodynamics [21 1. 

The goal of this work is to extract the shear viscosity coefficient r\ numerically from 
microscopic calculations using a standard setup motivated by the classical textbook picture 
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231 ] . In Fig. [TJwe introduce a particle system embedded between two plates. The two 



v z (x) 



-V 



wall 



FIG. 1: The classical definition of shear viscosity. Two plates moving in opposite directions with 
velocity ±i> wa ii. A flow gradient is established between the plates. The viscosity is proportional to 
the frictional force. 

plates move in opposite directions each with velocity i> wa u in z-direction. The moving walls 
are supplemented by two thermal reservoirs with ±t> wa u. In x-direction the system has an 
extension of size L. In y- and ^-direction the system is homogeneous and can be of infinite 
size. The mean free path of the particles should be very small compared to the system size, 
i.e. A m f p « L. On a sufficiently long time-scale a stationary velocity field v z (x) should be 
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established. In the non-relativistic limit the velocity field is linear. With the Navier-Stokes- 
ansatz the shear stress tensor TT tlu is proportional to the gradient of the velocity 

The proportionality factor is defined to be the shear viscosity coefficient rj. In Sec. [TT] 
we give basic definitions and information on the numerical model we use. In Sec. II III we 
demonstrate that Eq. ([1]) does not hold in general for the relativistic case, where the gradient 
is not necessarily linear and we discuss the shape of an ideal relativistic velocity gradient. 
Furthermore we will discuss the effect of viscosity and finite size effects on the velocity profile 
in Sec. IHH where an analytical formulation for the shape of the velocity profile is derived. We 
employ BAMPS to reproduce the velocity gradient as discussed in this chapter. In Sec. IIVI 
we compare the numerical results for the shear viscosity coefficient i] to an analytical value 
in order to confirm the applicability of our method. Finally we present the results on 
shear viscosity to entropy density ratio obtained from BAMPS with cross sections based on 
pertubative quantum chromo dynamics (pQCD) and compare them to existing calculations. 
We close with a summary. 

II. BASIC IDEA AND DEFINITIONS 

When systems are in stationary states, the first-order Navier-Stokes formulation of rela- 
tivistic viscous hydrodynamics can be used to calculate the shear viscosity 77, which is the 
proportionality factor between the shear tensor 7r^ u = and the velocity gradient V^it^: 



7T 



2r ] V < ^u u> , (2) 



where the projection 



\ (am + A M) - \ 



B aP (3) 



denotes the symmetric traceless part of the tensor B^ u . = — u^u v is the transverse 
projection operator and the metric is g^ v = diag(l, —1, —1, —1). 

Some definitions are in order. We use the Landau's definition of the hydrodynamic four- 
velocity 0: 

u^ 1 —^, (4) 
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where 

d 3 p 



is the energy-momentum tensor and the local energy density is defined as 

e = Ufl T^u u . (6) 

The shear tensor ir^ v is the difference of T^ v to its equilibrium value. For the geometry 
depicted in Fig. Q] = 7(1, 0, 0, v z ) with 7 = 1/ a/1 — v z . 

We will build up stationary states of particle systems via numerical simulations, which 
are realized by employing the microscopic transport model BAMPS, which solves the Boltz- 
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19]. In principle any 



mann equations for on-shell particles within a stochastic model 
microscopic transport model can be used for this purpose. 

Local values of tt^ and can be easily extracted from the numerical simulations by 
averaging over all particles contained in a bin of size Ax. However, to obtain the gradient 
of one has to take values from neighbouring local cells, which would cause additional 
numerical errors. To avoid such numerical problem we will first derive the analytical form 
of v z {x) for the given setup in Fig. [TJ Then we use this form and the numerically extracted 
7r^ to calculate the shear viscosity. 



III. VELOCITY, RAPIDITY AND FINITE SIZE EFFECT 

A. Analytical Derivation 

Instead of the hydrodynamic velocity v z (x) we address the position dependence of the 
rapidity y(x), which is defined by 

/ n 1 1 l + v z (x) ( 

y{x) = -ln- — . (7) 

2 1 - v z (x) 

Thus, v z (x) = tanhy(x). In the non-relativistic limit, where v z (x) is small, we have v z (x) ~ 
y(x). The advantage of y(x) is that it gets a shift by a Lorentz-boost e.g. with v z (xa) 

a m*a) [y( x )\ = y( x ) - v( x a) ■ (8) 



Demanding boost-invariance, i.e., A Vz ( XA )[y(x)} = y(x — xa), we obtain the solution y(x) = 
ax + b, where a und b are constant. Due to the boundary condition y(x = ±L/2) = ±y wa u, 
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y(x) is symmetric in x and thus, 6 = 0. If y(x) is continuous at the boundaries, we have 

y(x) = 2 -^x. (9) 

In the following we will convince ourselves from relativistic kinetic theory that Eq. ([9]) 
is only valid if the particle mean free path vanishes, or the distance L between two plates 
is infinitely long. For a non-vanishing mean free path and a finite distance L we will see 
discontinuities of y(x) at the boundaries. This is referred to as a finite size effect. 

We consider a general local observable A(x,t) with the definition 

A{x, t) = — V — J dT 1 F A (p 1 )f(p 1 ; x, t) , (10) 

where dT± = d 3 pi/(27r) 3 and n(x, t) = f dTif(pi;x,t) is the particle number density, p 
denotes the particle four-momentum. In our case n does not depend on position and time. 
In particular, for F A {pi) = np^/p\ we have the definition of particle four-flow A(x,t) = 
N^(x,t); for F^ipi) = \ ln[(p° +Pi)/(p? — Pi)] we obtain the rapidity A(x,t) = y(x,t) as 
given in Eq. ([7]), when using the Landau definition of the hydrodynamic four- velocity. For 
stationary states A(x, t) and the particle distribution function f(p; x, t) are constant in time. 

We define f(p; x, t) = f(p; x, t) / n(x, t), which is the probability density for the occurrence 
of a particle with momentum p around dT at (x, t). One obtains f(p;x,t) by summing 
probabilities for such events that a collision at (x', t') makes a particle having the momentum 
p and this particle travels to x at t without further collisions. It is mathematically expressed 
by 



f(pi,x,t) = 6(p lx ) / dx'wg ain (pi;x\t')wf ree (pi;x' \t';x,t) + 

J — oo 

POO 

O(-Plx) dx'w gain (pi]x' \t')Wfree(Pi;x' ',t']X,t) , (11) 



where Wg ain (pi; x', t') denotes the probability density that a particle with momentum pi is 
created via a collision at (x',f), and Wf ree (pi, x' ,t'; x,t) the probability that this particle 
travels from (x',t') to (x, t) without further collisions. Because f(p;x,t) is invariant under 
the transformation p — > —p, the two integrals in Eq. (TIT]) are equal. Thus, 

If 00 

f(pi,X,t) = - dx'Wgain{pi,X ,t )w free {p 1 \X ,t \X,t) . (12) 



2 . 

Our goal is to find the relation between A(x, t) and A(x', t'), which then can be used to solve 
A(x, t) analytically when the boundary conditions are given. 
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Using the standard definition of cross section for binary collisions of identical particles 

a22 = hJ §§ |MiM|2(27r)45(4)w +p * ~ Pi ~ P2) ' (13) 

where Myy-^xz is the matrix element and s = (pi + p 2 ) 2 = {p[ + P 2 ) 2 is the invariant mass, 
we have 

w gain ( Pl ; x', t')dx' = ^J dT' x dT' 2 f{p^ x', t')f(p' 2 ; x' , t')v rel ^dt' . (14) 

Vrei = s/(2pfp' 2 °) denotes the relative velocity for massless particles, dt' is the average 
time interval, during which a particle travels through dx'\ dt' = dx' < \p' x \/p'o > _1 and 
< \p' x \/p' >= 1/2 in thermal equilibrium. 

The probability Wf ree (pi, x', t'; x, t) is a product of w/ ree (pi; t"; x" + dx", t" + dt") over 
x" from x' to x: 

X X 

WfreeiPl] , t' ] X,t) = j [ W free (pi] x" , t" \ x" + dx" , t" + dt") = J| [1 - W loss (pi] x" , t")dx"} . 

x"=x' x"=x' 

(15) 

Wi oss {jp\] x", t") denotes the probability density that a particle with momentum p\ is destroyed 
via a collision at (x", t") and is expressed by 

WiossiPu x", t")dx" = J dT 2 f(p 2 ; x", t")v rel a 22 dt" , (16) 

where v re i = s/(2piP 2 ) and dt" = dx"(pf/p\ r )~ 1 . 

We now approximate wi oss (pi, x",t") to be the averaged one over p x : 

/I p x I 2dx" 
dT 2 f(p 2 ;x",t")(v re ia 22 )(^-)~ 1 dx" = 2n(v re ia 22 )dx" = - , 
Pi Amfp 

(17) 

where X m f p denotes the particle mean free path. This approximation applies for isotropic 
cross sections. In general, if the angular distribution is non-isotropic A m / P has to be replaced 
by an effective length scale, which is calculated as an average of the differential cross section. 
With Eq. (fT7|) we obtain the obvious expression 

/ 2dx"\ lx ~ x ' l/dx " ( 2\x-x'\\ 

w free (pi;x',t';x,t) = lim 1-t = exp . (18) 

dx"^0 \ A m f p J \ A m f p J 

Putting Eqs. flTJ), flU}, and (USD into Eq. (fTDJ) gives 

/°° 2\x-x'\ i r r ( j fT 

dx'e Wp - J dT^/ (pi; x>, t')f(p' 2 , x', t')v rel J dT 1 F A ( Pl )-^- . (19) 
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It is clear that replacing F^ipi) by Fa{pi) + Fa{j> 2 ) will leads to 2A(x,t). We now consider 
particular observables A(x,t) such that Fa is conserved in each collision, i.e., F^ipi) + 
Fa(p' 2 ) = F A (pi) + F A (p 2 ). We then have 

da 22 



/OO 2\x — x I 1 /* /* 

rfx'e~^7- y dV^dV'Jip^x'^fip^x'^FA^v^ J dTy^ 

OO 2\x — x'\ 1 /* /* 

rfx ' e -^7- J dT'^A^f^x'.t 1 ) J dT' 2 f(p' 2 ;x\t')(v rel a 22 ) 



2\x-x' 



dx'e x ™f? A(x' ,f) . (20) 



The same approximation is made as for wi oss in Eq. ffTT|) . Equation (TSUI) resembles the one 
derived in Ref. 



22j using "path integral method" in non-relativistic cases. 
We emphasize that Eq. (120]) holds only if the total Fa is conserved in collisions. For 
instance, the total particle velocity P1/-E1 + p 2 /E 2 is not conserved except in case the 
energy of all particles is same, whereas the total particle momentum rapidity is conserved. 
Therefore, the rapidity y(x) defined by Eq. ([7]) obeys Eq. (1201) . but the hydrodynamic velovity 
v z (x) does not. However, the total particle momentum rapidity is not conserved in 2 — > 3 
or 3 —7- 2 processes. In this case one has to take detailed balance into account and the sum 
of the total rapidity of a 2 — > 3 and its back reaction is conserved on average. If y(x) is 
conserved in collisional processes, it obeys Eq. f[2"0l . 

Equation (1201) represents a homogeneous first-order integral equation for A(x). It can 
easily be shown that the second derivative of A(x) vanishes, which leads to the solution 
A(x) = ax + b, where a and b are constant. We choose the boundary conditions to be 
A(x) = —y W aii for x < —L/2 and A(x) = y wa u for x > L/2 to reproduce the scenario 
indroduced in Sec. [B Since this scenario is symetric in x we have 6 = 0. To determine a we 
insert A{x) = ax into Eq. (|2"U|) and obtain a = 2y waU /(L + X m f p ). Finally the rapidity has 
the following form 

y( X ) = -^-X. (21) 

T A m f p 

We recognize the discontinuities of y(x) at the boundaries, which disappear only for vanishing 
mean free path A m / P — > or long distance L — > 00. Equation (12T|) is a new finding and 
accounts for finite size effects which must be taken into accout, if for numerical reasons 
X m fp/L cannot be made sufficiently small. 
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B. Numerical Confirmation 

In this subsection we will confirm our finding Eq. ( 12 ip by performing numerical transport 
calculations. We employ the parton cascade BAMPS. Details of numerical operations can 



be found in Refs. 



18 



191 ] . One important feature of BAMPS is that the model can simulate 
multiplication and annihilation processes such as the gluon bremsstrahlung process and its 
back reaction gg -h- ggg with full detailed balance. In order to verify the analytic findings 
we will first employ isotropic cross sections in BAMPS in the following. 

The numerical realization of the boundary conditions is as follows. Particles that reach 
the boundaries x = ±L/2 are removed, which simulates the particle absorption by the 
plates. Independent of the absorption, the plates emit particles, which pick up the velocities 
±v wa u of the plates. Here we treat the plates as thermal reservoirs of particles with the 
same temperature as those between the plates. The momentum distribution for emitting 
particles is proportional to the equilibrium Boltzmann distribution f wa ii{p) and the particle 
velocity p x /E: 

' 1X " n ~ -Uallip) (22) 



dtd?p E 
with 

M 

p M 11 

fwaii(p)=ge * , (23) 

where u^ all = 7 waM (l, 0, 0, v waU ), ^ wall = I/a/1 - v 2 wall , g = 16 is the degeneracy factor for 
gluons in SU(3), and T is the temperature. In the distribution (123"]) we neglect the quantum 
statistic factor for bosons and fermions. The rate of emissions can be calculated analytically 
(see App. [AJ and is 

^ = -^wall n wall , (24) 

where A wa u is the transverse area of the plates and n wa u is the particle density. In the xy- 
and xz-plane the boundary conditions are periodic. 

Particles between the two plates are initially distributed by the equilibrium form like 
Eq. (|2"3"]) with zero velocity. Figure |5] shows the buildup of the rapidity profile. We have 
chosen L = 2 fm, v wa ii = 0.5 (y wa ii = 0.55), and T = 0.4 GeV. Only binary collisions with 
a constant cross section are considered. The mean free path is set to be A m / P = 0.2 fm. 
Collision angles of the binary scatterings are distributed isotropically. 

The timescale for the buildup of the rapidity (or velocity) profile can be estimated as the 
mean diffusion time of particles traveling from one plate to another. For a gaussian diffusion 
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FIG. 2: Build-up with time of the rapidity profile. Results are obtained by averaging 500 events. 



process one has 



22 



< x 2 >= 2Dt. 



(25) 



where in the non-relativistic limit the diffusion constant D is the ratio of the shear viscosity 
7] to the mass density p. For relativistic case we replace p by the energy density e. As we 
will see in the next section, r\ w 1.2654nTA m / p = 0.42eA m / p (see Eq. fl29|) ). where e = 3nT 
is used. Thus, 

L 2 L 2 e L 2 
t = — = — « . (26) 



2D 2rj 0.82X mfp 

For our setup we find t ~ 24 fm/c, which is consistent with the numerical results shown in 

Fig. m 

Figure [3] shows the final rapidity profiles at sufficient long times. Calculations are per- 
formed for several mean free paths, in order to demonstrate the finite size effect. The 
differential cross sections are momentum-independent, i.e. isotropic. The lines present the 
analytical results given via Eq. (I2ip . while the symbols show the numerical values. One 
can see an excellent agreement, although approximations are made to obtain Eq. 02 II) . This 
indicates the validity of the approximations for using isotropic cross sections. 

On the contrary, when using the pQCD cross sections for gluons, which strongly depend 
on the invariant mass s one will have deviations from Eq. fl2T|) . The elastic and gluon 
bremsstrahlung process and its back reaction implemented in BAMPS are based on pQCD 
matrix elements given in Ref. |18[. Although the numerically extracted rapidity profile is 
different from the analytical form Eq. (|2"T|) . it is still linear in x. Replacing X m f p in Eq. 02 ip 
with an effective scale A e // one obtains the general formula. Using pQCD cross sections A e // 
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FIG. 3: Rapidity profiles for different mean free paths, A m f p = 0.02, 0.2, 2, oo fm. Constant cross 
sections and isotropic distribution of the collision angle are considered. The numerical results from 
BAMPS (symbols) are compared with the analytical ones (lines) given by Eq. (|2ip . 

has to be extracted numerically. Qualitatively, X e ff for pQCD interactions should be larger 
than X m f p , since pQCD-based processes prefer small angle scatterings and thus, are not as 
efficient for momentum transport as scatterings with isotropic angular distribution. 




IV. EXTRACTION OF SHEAR VISCOSITY 



In stationary states we can use the Navier-Stockes's formula Eq. (T2J) to calculate the shear 
viscosity rj. For the particular setup shown in Fig. [H Eq. fl2]) is simplified to 

7T" = -/J^M. (27) 

ax 

with 7 = 1/ a/1 — Vz(x). Using Eq. (1211 for v z (x) = tanhy(x) we obtain 

r, = -tt-VI - ■ (28) 

(Here X m f P replaced by A e //.) tt xz and v z (x) are extracted from BAMPS (the results of v z (x) 
are already shown in the previous section), whereas A e // is obtained by fitting y(x). 

Results for isotropic constant cross sections are presented in Sec. IIV Al Section IIVBI 
contains results for full pQCD interactions. 
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A. Elastic isotropic constant cross sections 



Elastic isotropic constant cross sections are meant that cross sections for elastic binary 
collisions are constant and the distribution of collision angle is isotropic. In__this case the 
shear viscosity of an ultrarelativistic Maxwell-Boltzmann gas is well known 



NS 



T 



1T ~ « i.zdm - = 1.2654 nT A mfp . (29) 
a 

serves as a benchmark to check the numerical methods applied to calculate 



Equation 
shear viscosity. 

Setups for numerical calculations are L = 2 fm, v wa u = 0.5, and T = 0.4 GeV. Results 
are averaged over N events = 2000 events. Figure H] shows the numerically extracted shear 
viscosity from Eq. ([28]) in each bin of size Ax = 0.2 fm. The mean value is in good agreement 



0.2 



BAMPS 
mean value 
analytic 



> 



0.1 



1 ■ 1 ■ 1 



x [fm] 

FIG. 4: Shear viscosity extracted from BAMPS and compared to the analytical result from Eq. (|29[) 
for A^p = 0.02 fm. 

with the analytical one from Eq. ( 129]) within the standard deviation, which decreases with 

1 / V 'events • 

Of couse this method for the extraction of shear viscosity can only be applied, when the 
particle system has relaxed to a stationary state. The relaxation time can be estimated 
according to Eq. ( [26]) . The relaxation time is inversely proportional to the mean free path 
and thus the extraction of shear viscosity for nearly perfect fluids with high N events is time- 
consuming. 

In Fig. [5] we show the mean shear viscosity with the standard deviation as a function of 
mean free path. The agreement with the analytical results (line) is perfect. This confirms the 
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FIG. 5: Shear viscosity as a function of mean free path. 

proposed method for extracting the shear viscosity for relativistic systems from numerical 
calculations. 



B. pQCD interactions 



In this subsection the results on the shear viscosity are presented for a system of gluons. 
For gluon interactions elastic gg — > gg and inelastic gg -H- ggg in leading-order pQCD based 



processes are included. For a more detailed discussion refere to Ref. 18|. 

Setups for this case are L = 40 fm, v wa u = 0.5, and T = 0.4 GeV. L has to be chosen 
appropriately as the mean free path increases with decreasing coupling constant a s . Running 
coupling is not implemented in the presented BAMPS calculations. 

The extracted mean values of the shear viscosity to entropy ratio rj/s are given in Table |T] 
and also shown in Fig. [7] with the resulting standard deviations of the simulations. The 
entropy density is taken by its equilibrium value s = An. 



TABLE I: r//s at various a s . 





0.01 


0.03 


0.1 
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0.3 
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0.6 




192.5 ± 23 


32.6 ± 3.49 


5.76 ± 0.63 


2.25 ± 0.2 


1.35 ± 0.14 


0.64 ± 0.064 


0.55 ± 0.06 




43.6 ± 5.2 


8.22 ± 0.66 


0.87 ± 0.09 


0.26 ± 0.03 


0.17 ± 0.01 


0.1 ± 0.01 


0.08 ± 0.01 
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FIG. 6: Shear viscosity to entropy density ratio at various a s 



0.01 



o.i 
_i 



BAMPS 
• Wesf 
Reininc 

Xu [PRL (2008) 100:172301 
El [PRC (2009) 79044914 



100 




0.01 



0.1 

pQCD coupling constant 



FIG. 7: Shear viscosity to entropy density ratio at various a s . Comparisons with other calculations 
using the same matrix elements for gluon interactions are made (more in text). 



The new results are compared to values by Xu et al. 20], El et al. 21] and Wesp et al. 



261 ] . Xu et al. identified the shear viscosity coeffient from the Navier-Stokes equation and 



used a gradient expansion in the Boltzmann equation. Calculating the second moment of the 
Boltzmann equation they obtained the shear viscosity coefficient in terms of the transport 
collision rate defined in Ref. 20j. El et al. derive the shear viscosity coefficient from the 



entropy principle, which also can be applied to derive the Israel- Stewart equations. For 
their derivation El et al. used Grad's approximation for the off-equilibrium distribution 
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function and obtained an expression for the shear viscosity similar to the one introduced by 

Xu et al., but with a sligthly different definition of transport collision rate. The results from 

I 

Wesp et al. [26] originate from equilibrium fluctuations. Here the Green-Kubo relations are 
employed to extract the shear viscosity. We see a good agreement with all three calculations 
for a s > 0.2. For a s < 0.2 our data is in exellent agreement with the data from El et 
al. as well es Wesp et al. We observe the 1/ (a 2 s log{l/ a s )) scaling behavior expected from 



Ref. 
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27]. 



Applicability of the methods by Xu et al. and El et al. crucially depends on the chosen 
parametrisations for the off-equilibrium distibution functions. In particular the momentum 
dependence of the off-equilibrium correction to the equilibrium distribution was chosen dif- 
ferently by these authors, which might explain the deviations between their results (Compare 



also the discussion in Ref. 



26|.) 



V. CONCLUSIONS AND OUTLOOK 



In this work we started with the classical picture of shear viscosity and shear flow. We 
demonstrated that the classical picture of a linear velocity field does not apply to relativistic 
systems. Rather, we found that velocity fields have a non-linear form, whereas the rapidity 
increases in fact linearly. 

We also derived an analytical expression for the rapidity and velocity profiles in systems 
where the mean-free path is non-zero. With an increasing mean-free path to system size ratio 
the slope of the rapidity profile decreases and finite size effects are not negligible anymore. 

We employed the numerical transport model BAMPS to create the velocity and rapidity 
profiles, compared the numerical results to our theoretical findings and observed an almost 
perfect agreement. The stationary gradient allows us to apply the relativistic Navier-Stokes 
equation to calculate the shear viscosity coefficient y. We found again a perfect agreement 
to the analytical value derived from kinetic theory |24j. The method proposed here to calcu- 
late the shear viscosity coefficient is thus perfectly suitable for other microscopic transport 
descriptions. 

Furthermore we have then used this setup to calculate the shear viscosity to entropy den- 
sity ratio in a numerical simulation with elastic and inelastic pQCD processes implemented 
in BAMPS for fixed coupling constant a s , which is varied from = 0.01 to 0.6. We compared 
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our results with previously published results [20|, |2l| and also with a very recent work based 



on the Kubo relation [26] and found a very good agreement. 
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Appendix A 

We calculate the number of particles AA" with a thermal Maxwell-Boltzmann distribution, 
passing through a wall of area A. The number of particles passing through a wall orthogonal 
to the ^-direction in a small timestep At is equal to the number of all paricles with distance 
Ax < —v T At from the wall: 



AiV = / -f^ / dydz [ dxf(p) = - [ -^Av x Atf(p) (Al) 

Jv x <0 l /7r J J A Jo<x<-v x At Jv x <0 l Z7r J 

AN f d 3 p 

-Av x f(p) 



At J Vx<0 (2tt 



, / d 3 p . , f d 3 p 



T 



d<p p t dp t / dyp t cos(4>)e t 

J TT JO J — OO 



,. . ... (A2) 

\ Z7V ) Jw J0 J -oo 

where w M = (cosh(/3), 0, 0, sinh(/3)). After the transformation of variables y — >■ y — (3 the 
dependence on the boost velocity drops out: 

AN gA f 2n ]A f 00 , f°° , wcchc,) 



AN gA f Z7T f°° f°° Pt c° S h 

aF = ^FI d( pj o p t dp t J_Jyp tC o^)e- r 



*2gA f 00 f°° o pt coghfa] 

dp 'L dme ~ T 



16 



2gA 




2T 3 gAT 3 nA 



cosh 3 (|/) (2tt) 2 ~ 4 



where n = gT 3 / tc 2 is the density in the local rest frame. 
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